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Research  Objectives 

(1)  Complete  the  development  of  a  measurement  system  designed  to  determine 
the  thermal  transport  properties  of  single  samples  of  either  zero,  one  or  two 
coating  layers  deposited  on  a  substrate. 

(2)  Measure  the  thermal  transport  properties  of  uncoatcd  metal  mirrors. 

(3)  Measure  the  thermal  transport  properties  of  single-layer  metal  coatings. 

(4)  Measure  the  thermal  transport  properties  of  two-layer  coating  stacks 
consisting  of  a  dielectric  layer  covered  by  a  metal  layer. 

Research  Status 

Objectives  1-3  were  successfully  completed.  Objective  4  was  not 
completed.  Appendix  I  contains  a  preprint  that  presents  a  comprehensive 
summary  of  the  research  results  obtained  under  the  present  AFOSR  grant. 

Publications 

The  principal  investigator  expects  to  submit  a  paper  describing  the 
currently  funded  research  to  Applied  Physics  Letters.  The  probable  title  is: 
’Thermal  Transport  Studies  of  Optical  Coatings,  Interfaces,  and  Surfaces.”  The 
anticipated  author  list  includes  the  principal  investigator  and  Gary  Wieinokly. 

Participating  Professionals 

This  research  was  carried  out  jointly  by  the  principal  investigator  and 
Gary  Wiemokly,  a  former  graduate  student  at  the  University  of  Southern 
California. 
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Oral  Presentations 


’’Thermal  Transport  Studies  of  Optical  Coatings,  Interfaces  and  Surfaces 
by  Thermal  Diffusion  Wave  Interferometry”,  Randall  T.  Swimm  and  Gary 
Wiemokly,  presented  at  the  1989  Boulder  Damage  Symposium,  held  Nov.  1-3, 
1989  in  Boulder,  Colorado.  Appendix  1  contains  a  preprint  of  the  paper 
submitted  to  the  conference  proceedings,  to  appear  under  the  title:  Laser  Induced 
Damage  in  Optical  Materials:  1989. 

Discoveries,  Inventions,  Patents 

No  discoveries,  inventions  or  patents  are  reported. 
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Appendix  1 


1989  Boulder  Damage  Symposium 


Thermal  Transport  Studies  of  Optical  Coatings,  Interfaces 
and  Surfaces  by  Thermal  Diffusion  Wave  Interferometry 

Randall  T.  Swimm  and  Gary  Wiemokly 
Center  for  Laser  Studies 
University  of  Southern  California 
Los  Angeles,  CA  90089-1112 

Thermal  transport  data  on  thin  metal  coatings  are  presented.  The  data  illustrate 
the  ability  of  the  measurement  technique  to  characterize  the  thermal  transport  properties 
of  single  samples.  The  method  offers  insight  on  coating  survivability. 

key  words:  adhesion;  coatings;  damage;  delamination;  laser  hardness;  laser-induced  dam¬ 
age;  nondestructive  evaluation;  nuclear  hardness;  optical  coatings;  thermal  conductivity; 
thermal  diffusivity;  thermal  properties;  thermal  transport;  thin  films. 

Introduction 

The  unifying  characteristic  of  thin-film  coating  studies  is  the  need  to  separate  the  surface  and 
interlace  properties  of  the  coating  from  the  bulk  properties  of  the  coating.  Frequently  this  separation 
is  accomplished  by  varying  a  characteristic  length  of  the  system.  Methods  of  doing  so  fall  into  at  least 
two  classes.  The  first  class  of  measurement  methods  is  based  on  varying  the  sample  thickness.  By 
plotting  the  measured  value  of  an  appropriately  chosen  property  of  the  samples  as  a  function  of  sam¬ 
ple  thickness,  the  sum  of  the  surface  and  interface  contributions  may  be  extracted  as  the  ordinate 
intercept  of  the  data  at  zero  thickness,  and  the  bulk  contribution  may  be  obtained  as  the  slope  or  the 
reciprocal  of  the  slope  of  the  data.  A  tacit  assumption  of  this  approach  is  that  the  sample  properties 
are  independent  of  coating  thickness.  An  additional  requirement  is  that  the  material  properties  be 
repeatable  among  samples.  Finally,  the  signal  to  noise  ratio  and  the  data  reproducibility  must  be  ade¬ 
quate.  In  practice,  it  may  be  difficult  to  satisfy  all  of  these  conditions,  especially  with  regard  to  mate¬ 
rial  properties  repeatability.  As  a  consequence,  data  are  occasionally  reported  which  show  nonphys¬ 
ical  results  such  as  negative  slopes  or  negative  intercepts.  Unless  the  substrate  polishing  and  coating 
deposition  processes  are  under  sufficient  control  to  yield  repeatable  coating  properties,  the  attempt  to 
obtain  meaningful  results  by  utilizing  measurement  techniques  based  on  multiple  samples  of  different 
coaling  thickness  is  futile. 

An  alternative  class  of  measurement  methods  is  based  on  varying  a  characteristic  length  im¬ 
posed  by  the  measurement  process,  and  obtaining  separate  surface,  bulk  and  interface  contributions 
for  a  single  coated  sample.  Such  techniques  may  include  those  based  on  varying  a  spot  size,  the  physi¬ 
cal  offset  between  pumped  and  probed  regions,  a  wavelength  or  a  diffusion  length.  The  present  work 
falls  into  this  latter  class  of  measurement  methods.  The  goal  of  this  paper  is  to  show  that  a  great  deal 
may  be  learned  about  the  thermal  transport  properties  of  a  single  sample. 
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Background 


The  measurement  and  analysis  techniques  have  been  described  in  earlier  proceedings*^. 
Briefly,  the  physics  of  the  measurement  is  as  follows.  Under  steady-state  conditions  of  sinusoidally 
modulated  surface  heating,  thermal  waves  generated  at  the  sample  surface  diffuse  through  the  sample 
and  partially  reflect  at  abrupt  discontinuities  of  the  thermal  effusivity  (equal  to  the  thermal  conduc¬ 
tivity  divided  by  the  square  root  of  the  thermal  diffusivity3).  Reflected  thermal  waves  that  return  to 
the  sample  surface  interfere  with  the  instantaneous  thermal  waves  being  generated  at  the  surface, 
causing  the  surface  temperature  response  to  be  out  of  phase  with  the  response  that  would  occur  in  the 
absence  of  reflected  thermal  waves.  In  the  present  work,  a  metal  coating  or  a  bare  metal  substrate  is 
used  to  ensure  the  required  condition  that  the  thermal  waves  be  generated  in  a  depth  that  is  small 
compared  to  the  local  thermal  diffusion  length  in  order  to  preserve  interference  contrast.  The  thermal 
waves  are  generated  by  heating  the  sample  using  a  sinusoidally-modulated  argon  ion  laser  beam.  The 
temperature  response  is  detected  using  an  unmodulated  helium  neon  laser  beam  reflected  from  the 
center  of  the  heated  region,  based  on  the  fact  that  the  metal  reflectivity  is  temperature  dependent. 
The  weak  intensity  modulation  induced  on  the  HeNe  laser  beam  is  measured  using  a  silicon  detector 
and  lock-in  detection.  The  system  electronic  phase  is  separately  calibrated  by  directing  the  pump 
beam  into  the  probe  detector.  The  measured  phase  of  the  temperature  response  constitutes  the  data. 
The  data  are  modelled  using  the  boundary  value  solution  for  the  presumed  sample  structure,  in  which 
the  material  properties  are  included  as  adjustable  parameters.  In  the  present  modelling,  the  piece- 
wise  constant  approximation  is  used,  and  the  thermal  contact  resistance  is  included.  A  key  considera¬ 
tion  is  the  uniqueness  of  the  solution.  This  can  be  inferred  only  by  individually  studying  a  sufficient 
range  of  different  sample  types  so  as  to  exclude  alternative  possibilities. 

The  present  measurement  system  is  an  upgraded  version  of  that  described  in  the  earlier  pro¬ 
ceedings*  '2.  A  Coherent,  Inc.  model  90-6  Innova  argon  ion  laser  is  used  as  the  pump  laser.  A  second 
probe  beam  line  has  been  added  to  allow  the  temperature  response  at  the  coating-to-substrate  inter¬ 
face  to  be  determined  for  samples  with  a  transparent  substrate.  Various  other  refinements  have  been 
added  to  avoid  systematic  errors.  These  include  improved  detector  shielding,  optical  baffling  to  mini¬ 
mize  light  scattered  into  the  detectors,  and  careful  optical  alignment  to  avoid  stray  reflections.  A  cus¬ 
tom-built  Lasermetrics,  Inc.  electro-optic  modulator  with  reduced  crystal  spacing  is  used  to  avoid 
thermal  lensing  in  the  index  matching  fluid  at  high  pump-beam  power  (3  W  at  514.5  nm).  Data  are 
recorded  using  an  automated  system. 

Recent  related  work  by  other  authors  is  listed  in  the  references'*"*  *. 

Results  and  Discussion 

Results  for  a  1  fim  rf  sputtered  Ni  coaling  deposited  on  fused  silica  are  presented  in  Fig.  I.  The 
sample  was  sinusoidally  heated  at  the  coating  surface.  The  phase  of  the  sinusoidal  temperature 
response  is  plotted  as  a  function  of  the  logarithm  of  the  modulation  frequency.  The  circles  show  the 
response  at  the  coating  surface.  The  squares  show  the  response  at  the  coating-to-substrate  interface. 
The  upper  curve  shows  the  results  of  a  zero-parameter  fit  (i.e.  no  adjustable  parameters),  based  on  a 
one-dimensional  heat  flow  model  assuming  the  handbook  values  for  thermal  conductivity  and  thermal 
dift usivity  of  Ni  in  bulk  format,  and  also  assuming  a  thermal  contact  resistance  equal  to  zero.  The 


lower  curve  shows  the  theoretical  response  at  the  interface  under  identical  assumptions.  Several 
conclusions  may  be  drawn  from  the  data  and  the  model  agreement. 


Fig.  1.  1  |im  RF  Sputtered  Ni  Coating  on  Fused  Silica 

First  consider  the  near-equality  of  the  surface  and  interface  data  at  low  modulation  frequencies. 
The  phase  of  the  temperature  response  at  the  surface  and  the  interface  is  expected  to  be  equal  at 
modulation  frequencies  for  which  the  coating  thickness  is  much  less  than  the  thermal  diffusion  length 
in  the  coating.  The  observed  near-equality  is  a  strong  indicator  of  whether  the  data  are  due  to  the 
temperature  modulation  of  the  coating  reflectance  as  assumed.  Because  of  the  small  values  typical  of 
the  thermal  coefficient  of  the  reflectance,  numerous  synchronous  background  signals  giving  rise  to 
systematic  errors  had  to  be  eliminated  before  the  observed  agreement  of  surface  and  interface 
response  was  attained. 

Attempts  to  fully  account  for  the  data  in  terms  of  a  one-dimensional  heat  flow  model  with  a 
non-zero  thermal  contact  resistance  failed.  However,  the  one-dimensional  model  is  known  to  fail  at 
frequencies  for  which  the  thermal  diffusion  length  is  not  small  compared  to  the  pump  spot  diameter. 
The  tentative  interpretation  of  the  data  is  that  the  Ni  coating  is  displaying  essentially  bulk-like  ther¬ 
mal  transport  behavior,  with  negligible  thermal  contact  resistance,  and  with  three-dimensional  heat 
flow  exhibited  at  low  modulation  frequencies.  In  order  to  verify  the  signature  of  three-dimensional 
heat  flow  in  the  simplest  possible  case,  a  solid  metal  mirror  was  studied. 

Results  for  a  solid  Mo  mirror,  using  a  1.90  mm  diameter  gaussian  pump  spot  are  considered 
next.  The  surface  of  the  mirror  was  visibly  oxidized.  Mo  is  known  to  form  sub-oxide  surface  layers 
several  hundred  angstroms  thick  12.  The  model  calculation  assumes  no  oxide  layer,  and  uses  hand¬ 
book  values  for  the  Mo  thermal  transport  parameters.  The  theoretical  derivation  follows  that  of 
Carslaw  and  jaeger^,  except  that  the  temporally  modulated  gaussian  pump  intensity  distribution  is 
used  for  the  heat  flux  instead  of  a  constant  flux  over  a  circular  region  as  treated  by  Carslaw  and 


Jaeger.  Consider  the  cyiindrically  symmetric  heat  diffusion  equation  for  a  homogeneous  isotropic 
medium  (assuming  small  temperature  excursions,  and  therefore  ignoring  temperature  dependence  of 
thermal  transport  coefficients),  for  the  region  z  2  0 


2 

a  t  a  T  i  aT 

y +  2  ~ - =  0  (1) 

ar2  az2  k  at 

where  T  equals  temperature,  ic  equals  thermal  diff usivity,  z  equals  depth,  r  equals  radial  distance  and 
t  equals  time.  For  a  solution  of  the  form  T  -  T  elut,  where  T  is  temperature  amplitude,  i  is  the  square 
root  of  minus  one  and  w  is  angular  frequency,  a  Fourier  transform  yields  the  steady-state  equation: 

2_  2_ 
a  T  a  T  loi 

ar7"  s?"7T'0'  ,2) 


For  a  gaussian  pump  beam,  Eq.  (2)  is  to  be  solved  subject  to  the  following  boundary  conditions: 


K  —  =  I0  e 

az  z-o  0 


where  K  equals  thermal  conductivity,  lo  equals  intensity  and  ro  equals  the  1/e2  radius  of  the  pump 
beam.  Applying  a  Hankel  transform  of  order  zero  to  Eq.  (2)  yields  the  equation: 

d2V«  2  iw 

— y  -  0  V0  _  —  V0  =  0  (4) 

dz^  K 


where 


r  T(r,z)  Jo  (or)  dr. 


and  0  is  the  transform  variable  conjugate  to  r.  Eq.  (4)  is  subject  to  the  following  transformed  bound¬ 
ary  conditions. 
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It  follows  by  inspection  that  the  following  is  a  solution  to  Eq.  (4),  subject  to  Eq.  (6)  and  vanishing  as 


1  (  -2(-J  ■  V" +!ir z 

v°  -  tt  I  re  Vro/  Jo(or )  dr  - —  .  - ... 
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and  therefore  by  inversion  that  the  following  is  a  solution  to  Eq.  (2),  subject  to  Eq.  (3)  and  vanishing  as 

z-*oo. 


Eq  (9)  expresses  the  temperature  at  a  general  point  in  a  semi-infinite  solid  with  gaussian  heating  at 
the  surface  In  the  experiments,  the  temperature  is  measured  at  the  center  of  the  heated  spot 
Therefore,  the  point  of  interest  has  coordinates  r=0,  z=0.  Consequently,  the  temperature  response 
simplifies  to  the  following  expression  (after  converting  the  integral  to  dimensionless  variables) *5; 


where  s  -  o^Jk/u),  and  b  -  (ro/2)Vw/2ic.  The  phase  ^  of  the  temperature  response  relative  to  the  sinu¬ 
soidal  pump-beam  modulation  may  be  obtained  as  the  arctangent  of  the  quotient  of  the  imaginary 
component  divided  by  the  real  component,  as  follows. 


Calculations  based  on  this  expression  (with  no  adjustable  parameters)  are  presented  in  Fig.  2. 

Reference  to  Fig.  2  shows  that  most  of  the  measured  phase  variation  is  accounted  for  by  the 
model.  On  the  basis  of  one-dimensional  calculations,  the  departure  of  the  data  from  the  model  is  con- 


sisteot  with  that  expected,  due  to  the  presence  of  the  surface  oiide  layer.  The  form  of  the  three- 
dimensional  numerical  calculations  also  supports  the  hypothesis  that  the  low-frequency  departure  of 
liie  measured  phase  from  the  one-dimensional  model  presented  in  Fig.  1  is  due  to  the  onset  of  three- 
dimensional  heat  flow. 


log  (f  (Hz]) 

Fig.  2.  Solid  Mo  Mirror 

In  an  effort  to  obtain  data  free  from  three-dimensional  heat  flow  effects,  a  metal  coating  of 
lower  thermal  diffusivity  was  studied.  Results  for  a  1  |im  rf  sputtered  Ti  coating  deposited  on  fused 
silica  are  presented  in  Fig  3  The  phase  of  the  sinusoidal  temperature  response  is  plotted  as  a  function 
of  the  logarithm  ol  the  modulation  frequency.  The  circles  show  the  response  at  the  coating  surface. 
The  squares  show  the  response  at  the  coating-to-substrate  interface.  The  upper  curve  shows  the 
results  of  a  zero-parameter  fit  (i.e.  no  adjustable  parameters),  based  on  a  one-dimensional  heat  flow 
model  assuming  the  handbook  values  for  thermal  conductivity  and  thermal  diffusivity  of  Ti  in  bulk 
format,  and  also  assuming  a  thermal  contact  resistance  equal  to  zero.  The  lower  curve  shows  the 
theoretical  response  at  the  interface  under  identical  assumptions. 

The  near  equality  of  the  of  the  surface  and  interface  data  at  low  modulation  frequencies  again 
confirms  the  absence  of  a  large  class  of  potential  systematic  errors.  Note  also  that  the  low  modulation 
frequency  data  are  in  agreement  with  the  one-dimensional  heal  flow  model  calculations.  The  gradual 
increase  of  the  phase  difference  between  the  temperature  response  at  the  surface  and  the  interface  is 
expected,  on  the  basis  of  thermal  transit-time  effects.  The  agreement  of  the  experimental  and  theo¬ 
retical  differences  between  surface  and  interface  response  per  se  is  fair.  Unfortunately,  the  data  show 
additional,  unexpected  structure  in  the  modulation-frequency  range  above  10  kHz.  The  authors  origi¬ 
nally  suspected  that  this  structure  may  have  been  due  to  the  presence  of  an  additional  stratified  layer 
with  an  abrupt  interface.  Possible  explanations  in  terms  of  a  thermally-thick  interface  layer,  a  coating 
surface  oxide  layer,  and  a  sub-surface  substrate  polishing  damage  layer  were  all  examined  via  mod¬ 
elling.  However,  additional  measurements  were  performed  on  a  number  of  samples,  including  a 


coaling  e  beam  deposited  onto  a  sapphire  substrate,  a  Ti  coating  dc  magnetron  sputtered  on  to  a  fused 
silica  substrate  that  had  been  argon-ion  beam  milled  after  polishing  to  remove  a  graded  layer  two  to 


Fig.  3.  1  fim  RF  Sputtered  Ti  Coating  on  Fused  Silica 

ten  micrometers  thick  16,  and  the  Mo  mirror  at  a  finer  grid  of  modulation  frequencies.  In  every  case, 
structure  was  observed  in  the  same  modulation  frequency  interval.  Since  the  unexplained  data 
structure  occurs  at  the  same  modulation  frequency  in  many  dissimilar  samples,  it  is  quite  likely  to  be 
the  result  of  a  systematic  error.  Further  work  to  identify  the  cause  is  under  way.  Notwithstanding  the 
presence  of  the  unexplained  structure,  the  tentative  interpretation  of  the  data  is  that  the  Ti  coating  is 
displaying  essentially  bulk-like  thermal  transport  behavior,  with  negligible  thermal  contact  resistance 
(much  less  than  10~2  K/(W/cm2)),  and  with  one-dimensional  heat  flow  exhibited  at  all  modulation 
frequencies  employed 

Conclusion 

Thermal  wave  interferometry  can  produce  a  wealth  of  thermal  transport  information  on  a  sin¬ 
gle  sample.  Onset  of  three-dimensional  heat  flow,  contributions  due  to  thermal  contact  resistance,  evi¬ 
dence  of  transit  time  effects,  and  the  presence  of  unanticipated  stratified  layers  are  some  of  the  phe¬ 
nomena  that  may  be  seen.  As  the  preceding  results  imply,  the  present  agreement  with  plausible  mod¬ 
els  based  on  handbook  transport  data  suggests  that  quantitative  modelling  using  adjustable  parame¬ 
ters  should  also  be  possible,  provided  that  the  physical  origin  of  the  observed  data  is  understood  The 
natural  consequence  of  the  sensitivity  to  the  detailed  transport  physics  that  is  shown  is  that  great  care 
must  be  taken  while  interpreting  the  data.  However,  the  revealed  detail  ensures  that  essential  effects 
are  not  overlooked. 
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Erratum 

A  systematic  error  that  entered  during  the  system  phase  calibration  was  discovered  in  the  data 
presented  in  (2|.  Consequently  the  data  and  data  analysis  in  that  paper  are  invalid. 

Addendum 

The  peak  in  the  data  shown  in  fig.  3  was  found  to  be  due  to  a  frequency-independent  2.5  degree 
error  in  the  orthogonality  of  the  in-phase  and  quadrature  channels  of  the  lock-in  detector.  The 
resultant  error  in  the  plotted  difference  data  (signal  phase  minus  calibration  phase)  shows  a  frequency 
dependence  according  to  the  actual  values  of  the  instrumental  phase  shift.  The  data  may  be  corrected 
by  a  simple  coordinate  transformation.  (8/21/90) 


